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Some rigorous results are presented for a first-order 
quantum phase transition between the dimerized state and 
\ Haldane-type state (i.e., a state similar to the ground state 
Q\ . of the one- dimensional spin-1 Heisenberg chain) in the spin- 
1/2 coupled chains with nearest- neighbour and next-nearest- 
neighbour Heisenberg interactions. Also presented are the 
£H . exact excited states in both phases. A partial phase diagram 
of the general spin-1/2 coupled chains is discussed. 

PACS numbers: 75.10.Jm, 75.50.Ee, 74.20Hi 
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Spin-1/2 coupled chains have attracted attention for 
the last few years [|l]-[| partly because of their relevance 
to such materials as (VO) 2 P20 7 @ and Sr 2 Cu 4 06 [||, 
and partly because of their intrinsic theoretical inter- 
est. For the isotropic coupled chains with only nearest- 
neighbour couplings, current consensus is that the system 
has a nonzero energy gap separating the ground and low- 
lying excited states, and the spin-spin correlation length 
of the ground state is finite (l]-|j|. No long-range order 
has been reported. 

It is well known that, for one-dimensional spin chain 
systems, the spontaneous dimerized phase (the perfect 
dimer state being defined as every two adjacent atoms 
forming a spin-singlet pair in the spin-1/2 or spin-1 
chains) || and Haldane-type phase (spin-1 chain) j?J are 
both possible. It is also well known that these two phases 
both exhibit nonzero gap and finite correlation length. 
The main difference between the two phases lies in the 
fact that the Haldane-type state is nondegenerate and 
translationally invariant with a hidden string order , 
while the spontaneous dimerized state of a single chain 
is doubly degenerate with translational symmetry bro- 
ken. For the system of spin-1/2 coupled chains, it is 
therefore natural to consider the possible existence of 
these two phases in despite of the dimerized phase be- 
ing non-spontaneous and to study corresponding phase 
transitions if the two phases do exist. In fact, spin-1/2 
coupled chains as a subject of research have a longer his- 
tory jlO 11 1 . The main interest there was to construct ef- 
fective one-dimensional chains with larger spin moments 
from spin-1/2 coupled chains. Many interesting results 
were produced [|lO|ll]]. In particular, it was found that 
spin-1/2 coupled chains can have the same ground state 
as that of the one-dimensional spin-1 Heisenberg chain. 



Because of their relevance to the experiments, the 
general spin-1/2 coupled chains are important models 
on their own right. In this article, we study spin-1/2 
coupled chains which clearly exhibit both the dimerized 
phase and Haldane-type phase when the coupling pa- 
rameters vary. Specifically, we consider spin-1/2 cou- 
pled chains with nearest-neighbour and next-nearest- 
neighbour Heisenberg interactions. The spin system con- 
sists of N pairs of S = 1/2 atoms, with interactions de- 
scribed by the following Hamiltonian, 



N 



r+1 



h = Y J i Js i - s 2 + J'csj-s; 



+ J"(Sl •S^ +1 + S' 2 -S^+ 1 )], 



S r 2 



(1) 



where the rung index r denotes each pair of spins with 
the spins on the top chain denoted as Si(r) and the 
spins on the bottom chain denoted as S2(j~); J is the 
coupling constant across each rung, J' is the nearest- 
neighbour coupling constant along each chain, and J" 
is the next-nearest-neighbour (diagonal) coupling con- 
stant. We mainly discuss the antiferromagnetic region 
(J, J', J" > 0). Periodic boundary conditions are taken 
in this article, namely, S^ +1 = S*, with fi = 1,2. The 
Hamiltonian of Eq. (1) at J = is the so-called compos- 
ite spin model JllJ . 

We first consider a single rung with a pair of spin-1/2 
atoms only. The rung has four states: the singlet state 
with zero total spin, namely |0) = -^=(| jj,) — | |f)) in the 
usual spin-1/2 notation; and the triplet states with total 
spin equal to 1, namely 1 1) = | TT>, |2) = ^(| TD + I IT)), 
and 1 3) = | ||) respectively. In a matrix representa- 
tion jl2],[l3| , these four states are denoted by four column 
matrices respectively, and the single spin operators (S* 
and with \i = 1,2) can then be written as 4 x 4 
matrices. We then introduce the so-called composite 
operator A nm (n,m = 0,1,2,3), which is a 4 x 4 ma- 
trix with a single nonzero entry, namely (n'\A nm \m') — 
Snn'dmm 1 - The single spin operators can then be writ- 
ten as linear summations of composite operators A nm . 
The meaning and bosonization of these composite op- 
erators have been discussed in detail in Briefly, 
Aqq (= 1/4 — Si • S2) is the singlet projection opera- 
tor, AQ n and A u q (n = 1,2,3) make transitions between 
the singlet and triplet states, and A nm (n,m — 1,2,3) 
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make transitions among the triplet states. A nm obey — X^r^ooi ^ s trivial, with no interaction between differ- 



pseudo-spin algebra, [A nm , A M ] = A n iS mk - A km 5 ln . 

In terms of these composite operators, A r nml with the 
index r denoting each rung in the coupled chains, the 
Hamiltonian of Eq. (1) can be rewritten as 
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+ (AI a + ^23) (A^ 1 + AU 1 } It is interest- 
ing to compare this Hamiltonian with that of the one- 
dimensional spin- 1/2 frustrated Heisenberg models (i.e., 
Eqs. (3.4)-(3.5) of Ref. Q). As we can see, by using 
the composite operators, the algebra of the spin- 1/2 cou- 
pled chains becomes quite transparent. It is trivial to 
note that when J' = and J" = 0, the ground state of 
the system is simply the perfect dimer state \D) which 
is given by the product of singlet states of all N rungs, 
namely \D) — Jl^Li More significantly, we notice 

that the Hamiltonian of Eq. (2) consists of three parts, 
the first part consisting only of the operators Aq which 
is nonzero only when acting on the singlet state of the 
r-th rung, the second part consisting only of the opera- 
tors H r>r+ i which makes transition between the singlet 
and triplet states, and finally the third part consisting 
only of the operators H' r r+1 which is nonzero only when 
acting on a state in which both the r-th and (r + l)-th 
rungs are in the triplet states. 

In fact, it is easy to prove that jJ^ r4 _ 1 is similar to the 
usual spin-1 Heisenberg interaction plj], namely 
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where we define operator P r by 

Pr = — ^33' Pr = V^C^ia + A 23 ), 
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with the usual SU(2) algebra, [P + ,P~] = 2P Z and 
[P Z ,P ± ]=±P ± . 

At J' = J", the second part of H disappears. The 
Hamiltonian is then reduced to 



H =J2(-A r oo+j'Pr-Pr+l), J" = J', 



(5) 



where, for convenience, we have set J' / J = j' and J = 1, 
and where we have ignored the constant term, N/4,. We 
notice that operators A^ and P r • P r +i commute with 
one another. We have therefore arrived at a Hamilto- 
nian consisting of two decoupled parts. The first part, 



ent rungs and with a gap value of 1 between the singlet 
and triplet states. The second part, J2 r P r -P r +i, is sim- 
ilar to the spin-1 Heisenberg chain with each rung in the 
coupled chains corresponding to each site in the spin-1 
chain (the only difference is appearance of a new type of 
excitation to be discussed later). The spin-1 Heisenberg 
model has been the focus of intensive investigations in 
the last decade and its physical properties are now well 
understood. The state-of-the-art calculations were re- 
cently performed by the density-matrix renormalization- 
group technique and extremely accurate results have 
been obtained. In particular, the long-range string order, 
<?(oo) = 0.374325096(2), is obtained by calculation of the 
string correlation function g(r) which is defined as 
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Since the two parts of H are decoupled and since all 
excited states of both parts have nonzero energy gaps 
(to be discussed later), its ground state is either that of 
the first part, or that of the second part, depending on 
the value of j'. The exact ground-state energy of H is 
therefore given by 



E, 
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-1 for f < j c , 
-f for j' > j c , 

Jc 



(7) 



where j c is the critical coupling constant, defined by 
j c = l/e — 0.713529353310(1) with eo the ground-state 
energy per site of the spin-1 Heisenberg chain |^] . Clearly, 
the phase transition at j' = j c is a first-order one because 
the first-order energy derivative with respect to j' is dis- 
continuous at the critical point. 

We next discuss excitation states in two phases. For 
the dimerized phase (j' < j c ), the excitations states 
can be easily obtained. For convenience, we define n 
singlet-to-triplet spin-flip (STSF) state as a configura- 
tion in which n rungs are in the triplet states, while the 
remaining (N — n) rungs are in the singlet state as in the 
dimer ground state. If each flip is separated from each 
another in a n-STSF configuration, the energy gap (i.e., 
total excitation energy minus total ground-state energy) 
of the n-STSF state is simply AE d (n) = n with 3™-fold 
degeneracy. However, if n flips form a contiguous clus- 
ter, the excitation gap can be much lower, with gap val- 
ues given by AEd(n) = n + E(n)j' , where E(n) are the 
eigenvalues of the open-ended n spin-1 Heisenberg chain. 
As well known, the ground state of the open-ended spin-1 
chain is singlet when n is even and triplet when n is odd; 
the energy difference between the singlet (triplet) ground 
state and next triplet (singlet) excited state for a given 
n decreases exponentially to zero in the large-n limit and 
the exact ground-state energy was obtained for up to 
n = 14 Q. For the one-flip state, E(l) = 0, therefore, 
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AEd(l) — 1 with triplet degeneracy. For the contigu- 
ous 2-STSF cluster, AE d (2) = 2(1 -j'), 2 -j', and 2+f 
for the singlet, triplet, and quintuplet states respectively, 
with a minimum gap of about 0.573 when j' — > j c . (For 
3- and 4-STSF clusters, the minimum gaps as j' — > j c 
are 0.895 and 0.685 respectively). If we consider a single 
STSF state as a (soliton) particle with spin momentum 
equal to one, it is obvious these particles attract to one 
another at low values of their total spin quantum num- 
bers. For large n (n > 14), E{n) = [— e (n — 1) + e'], 
where e' represents the residual boundary effect (its mag- 
nitude is expected to be much smaller than eo plj), the 
excitation gap of the n STSF cluster state is then given 
by 

AE d {n) = n+ [-e (n- 1) + e']f, n large (8) 

When j' — > j c , AE d (n) — > l+e'j c , a number independent 
of n. Since |e'| <§; 1, AE d (n) w 1 for large n. When n 
approaches iV, it is easy to see that the dimer ground 
state is unstable against the large n-STSF cluster state 
when j' > j c , as expected. 

The excitations in the Haldane-type phase have been 
discussed for the Hamiltonian of Eq. (1) at J = in 
Rcf. jy]. Here, we take advantage of the recent accu- 
rate results for the spin-1 Heisenberg chain Q, and fo- 
cus on the excitations of Eq. (5) for j' > j c . As we 
know, the low- lying excited state of the spin-1 chain is 
given by the one-magnon state with momentum q = n, 
and with the well-known Haldane gap A = 0.41050(2) 
||. Multi-magnon states have larger gaps HO]. We 
consider here a new type of excitations where the spin- 
1 chain is embedded with singlet rungs (llj. Again, 
for convenience, we define n triplet-to-singlet spin-flip 
(TSSF) configuration as a state in which there are n 
singlet rungs and the remaining (N — n) rungs are in 
their Haldane-ground-state configuration. We consider 
a single TSSF configuration, which corresponds to the 
ground state of the (N — 1) spin-1 Heisenberg chain 
with two open ends enclosing the singlet rung. Its to- 
tal energy can be written as Eh (I) = E(N — l)j' — 1, 
where, as before, E(N — 1) = — e (N — 2) + e'. Hence, 
AE h (l) = (2e + e')j' - 1. Since |e'| < 1, we obtain 
AEh(l) ~ 1 as j' — > j c , a value larger than 0.4105j c « 0.3 
of the one-magnon Haldane gap. Furthermore, from the 
numerical calculations of the open-ended spin-1 chains 
, we note that each of the two ends enclosing the singlet 
rung in the 1-TSSF state has an effective S — 1/2 spin 
with exponential decay of the local spin moment away 
from the ends, and with decay length £ w 6.03(1) (which 
is also the usual spin-spin correlation length). There- 
fore, the single TSSF configuration can be considered as 
particle-like (soliton) state with a size of about 14 lat- 
tice spacings (the core occupying two lattice spacings). 
We also notice that the degeneracy of the single TSSF 
state is four, while the one-magnon state of the spin-1 



chain is triplet. (It is also not difficult to obtain energy 
gaps of the 2-TSSF and other states and these details will 
appear elsewhere). Similar to the n-STSF cluster state 
in the dimer phase, the n-TSSF state has lowest energy 
when n flips are contiguous, with the gap 

AE h (n) = [(n + l)e + e']/-n. (9) 

We notice that when j' — ► j c , all clusters have the same 
energy gap of (1 + e'j c ) ~ 1, a number identical to that 
of the (large) n-STSF cluster state of the dimer phase at 
the critical point. When n N, Haldane-type phase is 
unstable against formation of the large n-TSSF cluster 
state when j' < j c , as expected. 

Finally, we discuss a partial phase diagram of the gen- 
eral Hamiltonian of Eq. (1). At J = J" = 0, the Hamil- 
tonian reduces to the two decoupled spin- 1/2 Heisenberg 
chains which is integrable by Bethe ansatz and which is 
critical with no gap. For J" = 0, we apply the (dimcr- 
ized) valence-bond spin- wave theory developed by us (ilf l 
via bosonization for the composite operators A T nm in Eq. 
(2), starting from the perfect dimer state \D) which is 
the exact ground state at J' = 0. Good results for 
the ground and spin-wave excited states are obtained for 
a small nonzero region of J' j J , but the theory breaks 
down at J'/J =1/2 and the dimerization order (similar 
to that of the one-dimensional spin-1/2 frustrated chain 
Hi) is nonzero in the region < J'/J < 1/2 @. This 
seems to suggest that the system of Eq. (1) without next- 
nearest- neighbour coupling (i.e., J" — 0) is dimerized for 
a nonzero region of J'/J and there is a critical value of 
J'/J at which the dimerized system is unstable against 
transformation to other phase(s). Combining with the 
rigorous results for J" = J' presented above, we propose 
a partial phase diagram as shown in Fig. 1. In Fig. (a), 
ABC is given by J" — J' with AB corresponds to the 
dimer phase and BC to the Haldane-type phase, and BD 
is only for guidance indicating the boundary between the 
dimerized phase and other possible phase(s). We denote 
the dimerized region enclosed by ABD as I. CE corre- 
sponding to the composite spin model is believed to show 
Haldane-type phase except the critical point at E [fTT|| . 
The dashed line between BD and E is quite arbitrary 
at this moment. It separates the Haldane-type phase re- 
gion (denoted as II) and any other possible phase region 
(denoted as III). 

There are three possible scenarios for III: (a) The 
valence-bond spin-wave theory is totally unreliable and 
the dimerized phase persists from I to III; (b) The 
valence-bond spin-wave theory is qualitatively correct 
and the whole region enclosed by BCED (i.e., II plus 
III) is Haldane-type phase, including the isotropic model 
with J = J' and J" = 0; (c) A new phase (such as 
spin liquid) appears in III which also has a nonzero gap 
but without any long-ranged order (dimerized or string 
order). We consider the first scenario most unlikely be- 
cause the (dimerized) valence-bond spin-wave theory is 



3 



expected to be at least qualitatively correct if the system 
has a nonzero dimerization order. Results from analytical 
approximations and numerical calculations seem to 
suggest the third scenario. However, since we have not 
seen any report on calculations of the string order, we 
cannot rule out the second scenario. Based on the rig- 
orous results along BC and numerical results along CE, 
we tend to believe the second scenario. In any case, it is 
interesting to apply powerful numerical techniques such 
as the density-matrix renormalization-group || to calcu- 
late the string order particularly in region III, using Eqs. 
(4) and (6). 

To conclude this article, we should also point out that 
there are other regions in Fig. 1 (a), where exact results 
can be obtained. For example, we show in Fig. 1 (b) the 
complete phase diagram for the Hamiltonian of Eq. (1) 
at J" = J', where GAB denotes the dimer phase region, 
BCF Haldane-type phase region, and FG ferromagnetic 
region. All of the three critical points, G, B and F, which 
separate the three phases, are first-order one. 
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